How much laser power can propagate through fusion plasma? 

Pavel M. Lushnikov 1,2,3 and Harvey A. Rose 3 
1 Landau Institute for Theoretical Physics, 
Kosygin St. 2, Moscow, 119334, Russia 
2 Department of Mathematics, University of Notre Dame, Indiana 46556, USA 
3 Theoretical Division, Los Alamos National Laboratory, 
MS-B213, Los Alamos, New Mexico, 87545, USA t 
(Dated: February 2, 2008) 

Abstract 

Propagation of intense laser beams is crucial for inertial confinement fusion, which requires 
precise beam control to achieve the compression and heating necessary to ignite the fusion reaction. 
The National Ignition Facility (NIF), where fusion will be attempted, is now under construction. 
Control of intense beam propagation may be ruined by laser beam self-focusing. We have identified 
the maximum laser beam power that can propagate through fusion plasma without significant 
self-focusing and have found excellent agreement with recent experimental data, and suggest a 
way to increase that maximum by appropriate choice of plasma composition with implication 
for NIF designs. Our theory also leads to the prediction of anti-correlation between beam spray 
and backscatter and suggests the indirect control of backscatter through manipulation of plasma 
ionization state or acoustic damping. 

PACS numbers: 42.65. Jx 52.38.Hb 
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INTRODUCTION 



Propagation of intense laser beams in plasma raises outstanding technological and scien 
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tific issues. These issues are closely tied with inertial confinement fusion (ICF) 
which requires precise beam control in order to maintain symmetry of spherical target implo- 
sion, and so achieve the compression and heating necessary to ignite the fusion reaction. ICF 
will be attempted at the National Ignition Facility (NIF). While most engineering features 
of NIF are now fixed, there are still crucial choices to be made |4J in target designs. Control 
of intense beam propagation is endangered by laser beam self-focusing, when a beam digs a 
cavity in plasma, trapping itself, leading to higher beam intensity, a deeper cavity, and so 
on. 

Self-focusing occurs when an intense laser beam propagates through a wide range of 



optica, m e dia L and na S oeen th e Subj ect of research fa m0r e th an forty yea,, since 
the advent of lasers (J. In laser fusion the intensity of laser beams is so large that self- 
focusing in plasma can cause disintegration of a laser beam into many small beams, leading 
to rapid change in beam angular divergence A8, called beam spray. Significant beam spray 
is absolutely unacceptable for attaining fusion which requires precise laser beam control 
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. It was commonly assumed that the main source of beam spray in fusion plasma is the 
self-focusing in local maxima of laser intensity (hot spots) which are randomly distributed 
throughout the plasma 4]. Hot spot self- focusing can be controlled by reducing beam 
correlation time, T c . However we show in this Article that the main limitation of maximum 
beam power, which can propagate in plasma without significant beam spray, is determined by 
collective instability which couples the beam to an ion acoustic wave. We call this instability 



collective forward stimulated Brillouin scatter (CFSBS) |7[ because it does not depend on the 
dynamics of isolated hot spots, but rather the intensity fluctuations as temporally smoothed 
(averaged) by ion inertia. We show below that this collective instability is consistent with 
the first experimental observation of the beam spray onset [^] while hot spot self-focusing is 
not. 
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BEAM COLLAPSE (CATASTROPHIC SELF-FOCUSING) 



There are two self-focusing mechanisms in plasma: ponderomotive and thermal. Histori- 
cally, pondero motive self- focusing was studied first. The ponderomotive mechanism results 
from averaging over fast electron oscillations in the laser electromagnetic field, at frequency 
ujq. Averaging induces an effective electrostatic potential proportional to the local laser 
intensity, which in turn adds to the usual fluid pressure term in hydrodynamical equations 
[jj. The thermal mechanism results from the transport of electron temperature fluctuations, 
ST e . 

Ponderomotive self-focusing in three dimensions (3D) is quite different than in two di- 
mensions (2D). (Here one dimension is the direction of propagation of laser beam with 
one/two transverse dimensions in 2D/3D, respectively). In 2D, self-focusing often results in 
propagation of optical pulses (called solitons [l(|) without change of their shape over large 
distances. In 3D, self- focusing often leads to dramatic intensity amplification with prop- 
agation distance. Indeed, self-focusing of light, as described by the nonlinear Schrodinger 
equation, results in formation of a point singularity after finite distance of light propagation 
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121 ]. A finite amount of optical power is drawn into this point, usually referred to as 



beam collapse. Near singularity, the nonlinear Schrodinger equation looses its applicability 
because of finite density depletion effects and instead of singularity, light scatters in a wide 
range of angles, causing loss of precise irradiation symmetry necessary for fusion. For appli- 
cation to fusion, only the 3D regime is relevant, and only this regime is considered in this 
Article. Note that in some regimes other, high frequency instabilities, such as stimulated 
Raman scatter can also arrest catastrophic collapse (see e.g. Ref. |13j) but they are not 
considered here. 

Beam collapse occurs if the laser beam power, P, exceeds a critical value 12], P c oc T e /n e . 



T e and n e are the electron temperature and density, respectively. For NIF parameters 
(n e « 10 21 /cm 3 , T e « 5keV, u Q « 5 x 10 15 sec -1 ) P c = 1.6 x 10 9 Watts. This power 
evaluation is based on Ref. fl2j| . in contrast to threshold given by Max jlJ ]. which is roughly 
half as large. The former may be dynamically realized (see Eq. (107) of Ref. Q]) from 
non-equilibrium initial conditions, appropriate to initiation by hot spots, while the latter 
is strictly an equilibrium property, and hence not useful for quantitative beam propagation 
prediction. 
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FIG. 1: Figure 1. Two dimensional slice of light intensity fluctuations inside plasma. Laser beam 
propagates from bottom of figure upward, (a) Distribution of fluctuations at scale much smaller 
than the beam diameter. Random fluctuations ("speckles"), are highly anisotropic, with correlation 
or speckle length along the beam propagation, "z", direction about 7F 2 Xq. (b) Beam spray regime 
of laser propagation. Beam disintegrates into many small beams, (c) Negligible beam spray regime. 
This regime is necessary for attaining fusion (from Ref. j^j with permission). Horizontal scale in 
(b) and (c) correspond to beam diameter. 



The energy required for inertial confinement fusion is so large that the power in each of 
NIF's 48 beam quads Q| exceeds P c by several orders of magnitude: the power of each NIF 
beam is approximately 8 x 10 12 Watts, or about 5 x 10 3 critical power. This difficulty is 
alleviated by the Random Phase Plate (RPP) which splits the laser beam into many 
(tens of thousands) small beams with random phases, which are then optically focused into 
plasma (see Figure 2 in Ref. [2]). As a result the total laser beam electric field amplitude, E, 
is well approximated in vacuum as an anisotropic random Gaussian field, with correlation 
length l c perpendicular to the beam propagation direction, much smaller than the parallel 
correlation length. The laser intensity, I oc \E\ 2 , forms a speckle field - a random in space 
distribution of intensity (see Figure la). 
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TIME-INDEPENDENT SELF-FOCUSING 



First consider the regime where laser beam time dependence is negligible. If the average 
intensity, (/), is small, then collapse events occur only in speckles (also referred to as hot 
spots) with / > (J), so that their power, P ~ I 2 1, exceeds P c . The width of these intense 
speckles, -FAo, is much smaller than the beam diameter and is determined by the laser optic 
system, where Ao is the laser wavelength in vacuum and F is the optic /—number (the 
ratio of the focal length of the lens divided by the lens diameter). We take l c = FXq/h. 
Since there is always finite probability of obtaining such collapsing speckles in the random 
Gaussian field model, the beam angular divergence, A8, increases with each collapse event. 
A9 in vacuum is given by A8 — 1/F, for F 2 ^> 1. If the probability of speckle collapse is 
small, then the beam will keep its initial form. But if laser power if so large that power of 
many hot spots exceeds P c then the beam will disintegrate into many small beams, leading 
to rapid change in A6, (beam spray). Figures lb and lc show examples of both regime of 
strong and negligible beam spray. 

An important measure of beam spray in this time independent regime is the fraction, 
Pscattered, of beam power, Pbeam, m speckles which self-focus as the beam propagates, es- 
timated as follows. NIF optic is approximately square, and hence a speckle area is F 2 \q, 
implying a critical intensity for speckle self-focusing, J c = P C /(FX ) 2 « 2 x 10 16 W/cm 2 . 

The a priori probability distribution of speckle intensities implies that the mean number 
M of speckles (local maxima) in volume V with intensities above value I is given by (see 
Eq. (21) of Ref. |l7|) 

jLV /2 _ a ( i_ x1I2 ~ 



M(I) ™ 



27F 4 A[j7r 



exp I ~ (7) I ' (1) 



.(I) J 

where (/) = Pbeam/ S is the average beam intensity, S is the beam cross section. Then 
M(/ c ) is the number of collapses per volume V and Pscattered = P C M(I C ) is the optical power 
scattered out of the main beam due to self-focusing. Therefore, rate of scattering is given 
by 



p-i dP Jdz P ° M &) (2 ) 

r beam UJr scattered / UZ — /Ay ' ^ ' 

For NIF parameters, Eqs. (0) and © give P h ~l m dP scattered /dz = 0.0002/cm for (J) = 
10 15 W/cm 2 and P b ~ l am dF \ C atteredi 'dz = 0.8/cm for (/) = 2 x 10 15 W/cm 2 . If Max's lower value 
of P c were used, order unity of the total beam power would have been predicted to scatter 



over a typical NIF plasma length of 1cm, even at the lower intensity since dP scattere( i/dz is 
exponentially sensitive (see Eq. ((H)) to the parameter a, with dP sca ttered/ dz oc exp(— a), and 
a = Ij (I). For NIF parameters, a 1. 



TIME-DEPENDENT SELF-FOCUSING 

Clearly beam spray due to speckle self-focusing could be a problem at the higher intensity. 
This is alleviated by temporal beam smoothing techniques jlfiillfll which induce finite speckle 
coherence time, T c : the intensity distribution of light intensity inside plasma is given by a 
speckle field at each moment of time as in Figure la but location of hot spots changes 
in a random manner with a typical time T c . Such techniques are used in contemporary 



experiments 



and in future experiments at NIF. 



Inertia weakens the plasma density response: if T c is less than the duration of a particular 



self-focusing event, ~ FX /c s ^/ P/ P c , (this estimate is accurate for P/P c ^ 2.5, see Ref. 
jl5^ ) then this self- focusing event will be suppressed. This suppression effect is significant 
if T c < FX /c s , i.e., T c must be smaller than the time it takes for a sound wave to cross a 
speckle width (~ 4ps for NIF parameters). Here c s is the ion-acoustic wave speed. (This is 
in contrast to the case of almost instantaneous response of optical Kerr nonlinearity which is 
typical for solids [5j). As T c decreases, a smaller fraction of the beam power participates in 
collapse events, controlled by the parameter a{l c / c s T c Y, instead of a, for time independent 
self-focusing. This has led to the common assumption 4| that if the total power participating 
in independent collapse events is made arbitrarily small by reducing T c , then beam spray 
could be reduced to any desired level. 

However, we have found |7[ that even for very small T c , self-focusing can lead to strong 
beam spray. Now, self-focusing results from a collective instability, CFSBS, which couples 
the beam to ion acoustic waves that propagate transversely to the direction of laser beam 
propagation. As Lincreases, the well-known dispersion relation of forward stimulated Bril- 
louin scattering jjjj] is recovered for coherent laser beam. We predict that this instability 
is not a sensitive function T c for c s T c < FXq. Recent experiments at the Omega laser fa- 
cility [8] are in excellent agreement with that prediction: It was found that reducing T c 
from 3.4ps (for which c s T c w FXq) to 1.7ps did not cause a further reduction of beam 
spray at (/) = 5 x 10 14 W/cm 2 . Note that dominant seed for CFSBS is not thermal but 
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time-dependent plasma density fluctuations caused by fluctuating speckles. 



THERMAL SELF-FOCUSING 



n 



Quantitative comparison with this data requires extension of our earlier work [7] to allow 
transport of fluctuations, ST e , in electron temperature. In that case the second mechanism of 
self-focusing - thermal self-focusing comes into play. Propagation of laser beam is described 
by paraxial equation for the electric field spatiotemporal envelope, E, 



which is coupled to linearized hydrodynamic equation for the relative density fluctuation, 
p = 5n e /n e , as it propagates acoustically with acoustic speed c s : 



where ST e is the fluctuation of electron temperature, ko = 2tt/\q, I = \E\ 2 is the light 
intensity, v is an integral operator whose Fourier transform in x and y is Ui a kc s , where vi a 
is the ion acoustic wave amplitude damping rate normalized to the ion acoustic frequency. 
x and y are transverse directions to beam propagation direction z. E is in thermal units 
defined so that in equilibrium, with uniform E, the standard p = exp(— /) — 1 is recovered. 
n c = m e ujQ/ '47re 2 is the critical electron density, m e is the electron mass and e is the electron 
charge. The relative electron temperature fluctuation, Sr.fi. is response for therma, 
self- focusing and was omitted in our previous work p|. 



We make the ansatz that the Fourier transform of electron temperature fluctuation, 
ST e (k)/T e satisfies, 



which is a reduced version of Epperlein's model [25]. Here the right-hand-side (r.h.s.) de- 
termines plasma heating by the inverse bremsstrahlung, I(k) is the Fourier transform of /, 
so that intensity fluctuations are a source of 8T e jjj. The inverse bremsstrahlung relaxation 
time, Tib, is given by, 
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Also 

_ [i + (3oa e )V3] 



and is an empirical factor |23|, 4> = (4.2 + Z*)/(0.24 + Z*), Z* = E™^ s VXX^i is the 

i i 

effective plasma ionization number, rij and Zj are the number density and the ionization 
number (number of ionized electrons per atom) of i-th ion species of plasma, respectively. 
A e is related to the standard e — i mean free path, X ei , by A e = (X ei /3)(2Z* /ii(p) l/2 . The 
basic ion acoustic wave parameters, u ia and c s , are regarded as given by kinetic theory 
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28, 



29| which, e.g., takes into account the effect of compressional heating on sound wave 
propagation. For comparison with experiment in this paper, however, collisionless theory is 
used for evaluation of acoustic wave parameters. 

Eq. (|SJ) implies that thermal conductivity is determined by 

- 3 Ue - ksH (Q) 

K ~ 2 7-aJfe 2 ~ l + (30A;A e ) 4 / 3 ' 1 j 

where ksH is the classical Spitzer-Harm jsjj] thermal conductivity coefficient in plasma. 
Since l c is not small compared to the electron ion mean free path, X ei , thermal transport 
becomes nonlocal, and ksH is effectively reduced, as given by Eq. ©, when applied to a 
fluctuation at speckle wavenumbers, k = 0(l/l c ). This reduction of ksh is substantial for 
experiment of Ref. implying much larger 5T e than classical transport 25]. Importance of 
ihe thermal contribution to self-focusing at the speckle scale was first realized by Epperlein 

in ~i 

22, 25 , on the basis of Fokker-Planck simulations, and later analytically derived 24 and 

n nn \' 

verified experimentally j26|. It was recently realized [27, 28] that Epperlein's result |2J 
is correct provided the acoustic frequency c s /l c is smaller than the electron- ion collision 
frequency v e / X ei . 

To solve Eqs. (jSJ), © and JHJ) we need to determine boundary conditions on E. We 
assume, absent plasma, that in the optic far field the Fourier spectrum of E is top-hat with 
square shape: 

|J5(k)| = const for < k m and \k y \ < k m ; \E(k)\ = 0, otherwise, (9) 

where k m = l~ x . Thus our boundary conditions correspond to square top hat. The superpo- 
sition of all these Fourier modes propagating in uniform density plasma we refer to as Eq, 
the solution of Eq. (j3J) with p = 0. We assume temporal beam smoothing which means that 
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Fourier modes E(k) with different k are uncorrelated and the modes with the same k are 
correlated with short correlation time T c < l c /c s . 

For NIF designs, Z* is highly variable depending on details of plasma composition. Laser 
beam may pass through, e.g., He, Be, CH , S1O2 and Au plasma, allowing a wide range of 
Z*. When Z* is small, thermal effects are small, and our previous ponderomotive theory 
0] applies. In this case, the linear stage of the collective instability depends only on one 
parameter - dimensionless intensity 

AF 2 n P 1 



'0 



v ia n n au, 

12 



/o oc . (10) 



10 is the spatial average of \E\ . Note that the standard figure of merit for self-focusing, 

1 1 a, is smaller by the factor v ia (see Ref. jj]). 



COLLECTIVE FORWARD STIMULATED BRILLOUIN SCATTER AND TRAN- 
SITION TO BEAM SPRAY REGIME 

For small T c , one might expect p ~ and that the laser beam would propagate with 
E = Eq. However, linearization of Eqs. (JHJ), (@J) and © about this state shows that this 
propagation is unstable. Following ideas of Ref. j| and setting p = 5pe Xz exp [z(k • x — u;i)] , 
E = Eq + 5Ee Xz exp [i(k • x — ut)\, we obtain the following dispersion relation, at acoustic 
resonance uo = kc s , assuming k parallel to either the x or y directions: 



1 g(k\ t 



I (id 

dp 



1 ikc s Tii) 

where the plasma density response function 51 /5p is given by 

5J_ = n 1 tfI^ k\-2k m + kf + Akl\ 2 

5p n c 4kk m n k 2 (2k m + k) 2 + Ak 2 X 2 ' 1 j 

In general case of arbitrary direction of k the dispersion relation is much more bulky and 
not given here because it gives essentially the same result. 
Note that for k m -»• (F 2 > 1) Eq. (JEJ) reduces to 

6I_ _ n e 2k 2 kll 

5p n c 4k 2 \ 2 + k* [ ' 

which means that Eq. (|TTj) . absent thermal effects (i.e. for ST P = in Eq. ((jlj)), reduces to 
the paraxial limit of the standard FSBS dispersion relation [2 



Absent thermal effects we regain the pondermotive case considered in Ref. (tJ except 
that in this Article square top hat boundary conditions Q are used compared with circular 
top hat boundary conditions used in We find however that both circular and square top 
hat boundary conditions give similar results. 

Positive value of Re(X) corresponds to convective instability so that the fluctuations 
of beam intensity grow as exp [i?e(A),z] with distance. When A is non-dimensionalized, 
A = l 2 k X, it only depends on 

g(k\ e ] 



1 + 



1 ikc s Tib 



Iq. (14) 



Here k is determined from the condition that Re(X) is maximum and Iq is given by (jlOj) . 
According to our theory of CFSBS, A -1 , should be compared with the basic correlation length 
in z direction, known as the speckle length, l spec kie ~ 7F 2 \q. The value A = 0.1, at which 
A ~ hpeckiei marks regime transition. In the first, weak regime, with A <C 0.1, there is little 
gain over a speckle length. It follows that only small changes in correlations develop over 
a speckle length, in particular, there is little change in A6. Changes over different speckles 
are uncorrelated, leading to a quasi-equilibrium (see Figure 3 of Ref. [7|). As A crosses the 
value 0.1 (corresponding to / ~ 2 in ponderomotive case), a second, non-equilibrium regime, 
is entered, and beam properties change rapidly with z. In particular, A8 changes rapidly, 
i.e., there is beam spray. This is shown in Figure 2, where normalized beam spray rate 
is shown. Note that absent instability one expects beam spray rate d(8 2 )/dz ~ Iq. So in 
Figure 2 we normalized beam spray rate to Iq (see Ref. j3| for more discussion). Compared 
with Figure 6 of Ref. Q|, there has been an important change of independent variable, from 
/ to A, which allows a unified presentation of both ponderomotive and thermal cases. 

Thus analysis of A results in second and main conclusion of our CFSBS theory: predic- 
tion of the onset of beam spray, and hence a prediction of fundamental limit on power 
propagation. Here we present comparison of this prediction with M, the first experi- 



mental measurement of beam spray onset (see Figure 3). From [jc], [n| we find that 
0.14 < n e /n c < 0.25. T e ~ 2keV, F = 6.7, uj « 3.6 x lO^sec" 1 , and Z* = 6.4 at upper 
range of densities. For a nominal electron density of n e = 0.2n c , the 0.1 contour (color 
online) of A is shown in Figure 4a, implying / w 0.65 at regime transition, was observed 0, 
corresponds to / ~ 1.05, with Landau damping v ia = 0.06 for the plasma composition at 
this density. The major uncertainty in comparing this data with theory is due to significant 
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(a) 




(b) 
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FIG. 2: Figure 2. Dependence of beam spray rate obtained from simulations as a function of 
(a) dimensionless growth rate A and (b) dimensionless intensity /. Red curves correspond to 
ponderomotive self-focusing {Z* = 0) and black curves (Z* = 15) correspond to case there both 
ponderomotive and thermal self-focusing are essential. Both red and black curves collapse to the 
single curve in (a) which indicated that A is a much better parameter for onset of beam spray, 
compared with /. 



time dependence of T e /Tj during experiments as well as plasma density inhomogeneity, e.g., 
if n e = 0.14n c (which corresponds to plasma density plateau in Figure 3 of Ref. jsj]) with 
other parameters the same, then theory predicts / ~ 0.73 and experiments give / ~ 0.82. In 
contrast, prediction based on speckle collapses, gives that even at the maximum density of 
n e /n c = 0.25, P^ am dP sca tteredi 'dz = 0.23cm" 1 , the scattered power fraction, F Mattered / 'Pbeam, 
is only 0.5% after 200/zm of propagation through the high-density region of the plasma. This 
is much less then the observed |a] 10%. Therefore, beam spray due to CFSBS is consistent 
with the data while beam spray due to speckle collapse is not. 



IMPLICATION FOR BACKSCATTERING 



Recent experiments at the Atomic Weapon Establishment in the UK have demonstrated 



reduction of both stimulated Brillouin and Raman backscatter 



32| by the addition of small 



amounts of high ionization state dopants to a low ionization state plasma, e.g., a 1% dopant 
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(a) 



(b) 




Figure 3 

FIG. 3: Figure 3. Experimental images of cross section of time averaged laser beam intensity after 
propagation through plasma, (a) Onset of beam spray regime at 5 x 10 14 W/cm 2 . (b) Negligible 
beam spray regime achieved by lowering intensity. Dashed circles correspond to F = 6.7 beam 
width for propagation in vacuum. Reproduced from J] with permission. 

reduced backscatter by more than an order of magnitude. Combination of these experimen- 
tal facts with our prediction that dopant may cause transition to beam spray regime suggests 
that one should expect anti-correlation between beam spray and backscatter. If this anti- 
correlation is confirmed experimentally then we propose the following mechanism: beam 
spray decreases speckle length (correlation length) with beam propagation and backscatter 
is suppressed by reduction of laser beam correlation length. The latter has been estab- 
lished through simulation [13|, experiment and one dimensional analytic theory 34]. 
In other words, control of backscatter is achieved indirectly through control of CFSBS. We 
are unaware of any other explanation of this backscatter reduction by the addition of small 
amounts of high Z dopant. 

Clearly, to maintain control of forward beam propagation, beam spray must not be strong. 
If plasma paramters are conducive to backscatter as in the Atomic Weapon Establishment 



experiment 



then by altering the plasma state so as to be above, but close to, the beam 



spray regime transition, allowing moderate beam spray might lead to optimum control of 
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FIG. 4: Figure 4. (a) Solid curve separates predicted beam spray regime, A > 0.1 (from green to red 
colours), from negligible beam spray regime A < 0.1, (from blue to purple colours). Different colours 
denote values of A, with red corresponding to the value 0.3 and above. Magenta square denotes 
experimentally measured (Ref. ^|) beam spray onset, assuming Vi a = 0.06 and black circle is the 
theoretical prediction for z/j a = 0.06. (b) Predicted onset of beam spray regime (i.e. for A = 0.1) as 
a function of Z* and Vi a for NIF plasma with T e ~ 5keV, F = 8, n e /n c = 0.1, u>o ~ 5.4 x 10 15 sec _1 . 
Colours show laser intensity, in units of 10 15 W/cm 2 . Intensity is at maximum for small Z* and 
large i/j a . We assume Z* > 4 to make sure that condition c s /l c < v e /X e i is true. 



beam propagation and backscatter. This suggests operating above but, e.g., close to the 
solid curve of figure 4a which marks the transition regime of CFSBS. 



CONCLUSION 



In conclusion, transition to the beam spray regime was recognized as a collective phe- 
nomenon. Our theory is in excellent agreement with experiment: the transition laser inten- 
sity and its insensitivity to changes in correlation time were predicted. We found that the 
growth rate of CFSBS depends on four dimensionless parameters: the scaled laser intensity 
J (see Eq. ()ll)jl ). scaled electron-ion mean free path A ei /FA , effective ionization number 
Z* , and c s /v e . The first three of these can be manipulated experimentally. So our theory 
permits predictions for beam control at NIF that may be implemented since thermal self- 
focusing can be manipulated experimentally through control of CFSBS in two ways. First, 
by changing Z* through change of plasma composition. For example, addition of 1% of 
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Xenon (high Z dopant) to low Z plasma (50% of He and 50% of H) would increase Z* from 
1.7 to 15.5 without significant change in Ui a . Second, beam control can be implemented by 
adding low Z dopant to a high Z plasma, e.g., adding He to Si0 2 , in order to increase v ia 
at almost constant Z*. Figure 4b shows dependence of laser intensity (indicated by colors) 
at predicted onset of beam spray regime on Z* and v ia for NIF parameters. It is seen that 
maximal allowable intensity occurs for small Z* and large v^. We propose Figure 4b as 
direct guide for choice of NIF designs to attain maximum power of laser beam, which may 
propagate without significant beam spray. 

Observation of anti-correlation between beam spray and backscatter, through the addtion 
of small amounts of high Z dopant, would mean addtional confirmation of our theory. We 
predict that control of backscatter is achieved indirectly th roug h control of CFSBS, e.g. by 
changing plasma ionization state and/or acoustic damping 35j . 
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